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ABSTRACT 

We report the discovery of 73 double-mode RR Lyrae (RRd) stars in fields 
near the bar of the LMC. The stars are detected among the MACHO database 
of short-period variables that currently contains about 7900 RR Lyrae stars. 
Fundamental periods (Pq) for these stars are found in the range 0.46-0.55 days 
and first overtone-to-fundamental period ratios are found to be in the range 
0.742 < Pi/Pq < 0.748. A significant fraction of our current sample have period 
ratios smaller than any previously discovered RRd variables. We present mean 
magnitudes, colors, and lightcurve properties for all LMC RRd stars detected to 
date. The range in period ratios is unexpectedly large. 

We present a determination of absolute magnitudes for these stars based 
primarily on pulsation theory and the assumption that all observed stars are 
at the fundamental blue edge (FBE) of the instability strip. Comparison of 
the calibrated MACHO V and Kkc photometry with these derived absolute 
magnitudes yields an absorption-corrected distance modulus to the LMC of 
18.57 ± 0.19 mag which is in good agreement with that found (18.5) through 
comparison of galactic and LMC Cepheids. 

Adopting this luminosity calibration, we derive an increase in the distance 
modulus, and thus a reduction in the age found via isochrone fitting for M15 of 
about 33% and discuss the implications for cosmology. 

Subject headings: stars: RR Lyrae : LMC 
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1. Introduction 

Multimode RR Lyrae, also known as RRd type, are especially interesting tests of 
pulsation and evolutionary theory because precise period ratios allow an estimate of the 
mass of each star. This paper reports properties of a total of 73 such stars in the Large 
Magellanic Cloud (LMC), nearly twice as large as the number of all known RRd stars prior 
to this work. 

The MACHO Project LMC RR Lyrae database is described in Alcock et al. (1996) 
where accumulated V and R photometry for approximately 7900 RR Lyrae stars in 22 
intensively-observed fields has been analysed to determine the distributions of period and 
amplitude. In that paper, a random selection of 500 lightcurves was examined and two 
multimode stars were identified. The ratios of first overtone-to-fundamental period, Pi/Pq, 
of those stars were found to be lower than any previously discovered RRd stars. 

This paper reports the discovery of a further 71 RRd stars. We first discuss the 
observational properties of this new sample. Next we estimate absolute magnitudes using 
pulsation theory. Finally, we note implications of the derived luminosities for cosmology. 

2. Results 
2.1. Photometry 

The photometric reduction procedures for the MACHO images are described in Alcock 
et al. (1996). The MACHO instrumental magnitudes have been placed on the standard 
photometric system of Kron-Cousins V and R through comparison with Landolt (1992) 
standard stars that were transferred to the 22 high-priority MACHO fields. While gross 
zero-point differences between fields have been removed, small offsets within each field exist 
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as an artifact of the MACHO photometry reduction procedure. We conservatively adopt an 
error in V or R of ±0.10 mag. The transformation between MACHO instrumental color and 
Kron-Cousins (V-R) is less affected by the internal zero-point offsets, therefore, the color 
for each star in Table 1 has an estimated uncertainty of ±0.05 mag. We arrive at these 
estimates of the photometric uncertainties through numerous consistency checks within the 
MACHO database and by comparisons with published photometry of LMC stars. 

The lightcurves used in this analysis represent all data acquired for these stars during 
an interval of approximately 1400 days and typically represent 400-700 points per color. In 
Figure 1 we show the mean apparent V and (V-R) for LMC double-mode RR Lyrae and 
illustrate their position relative to other LMC field stars. The lightcurves for these stars are 
available on the MACHO Project home pages at the URLs [http : / /wwwmacho . anu . edu . au/ 



and http : //wwwmacho .mcmaster . ca/ . 



2.2. Selection Criteria 

Double-mode RR Lyrae candidates were selected from our color-magnitude diagram- 
based RR Lyrae sample by two techniques. The first, described in Alcock et al. (1996), 
relied on visual examination of 500 lightcurves using the IRAF PDM task. More recently, 
candidate RRd stars were selected based on their measure-of-fit to a singly-periodic 
lightcurve, their period and lightcurve amplitude. 

The strong period dependence of our selection criteria deserves mention. For 
photometric periods between 0.35 and 0.40 days - corresponding to fundamental periods 
between 0.46 and 0.50 days - our selection criteria identify double-mode RR Lyrae very 
efficiently. Shortward and longward of this range, our yield is reduced by the increase 

^^Developed at the National Optical Astronomical Observatories 
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in false positives with period estimates of one-third and one-half of a day, respectively. 
The measure-of-fit statistic is less effective at long periods due to the presence of Blazhko 
variation in some RRab pulsators. Therefore, we caution the reader that our current sample 
is not complete and that the distribution in period (and therefore mass) is unrepresentative. 

2.3. Observed Properties of the Sample 

We summarize our discovery of LMC RRd stars in Table 1. This table lists, from 
left-to-right, the equinox J2000.0 right ascension and dechnation, the internal database 
identifier, the number of measurements and mean magnitude in V and Rrc, respectively, 
the estimate of the fundamental period Pq (in days) and the period ratio, Pi/Pq, the 
semi-amplitude of the Fourier component for the frequency corresponding to Pi (in mags), 
and the ratio of the amplitudes for the first Fourier components of the first-overtone and 
fundamental modes. The uncertainties in the period and period ratio are expressed in 
parentheses in units of the least significant digits. Table 1 contains 75 entries. The pairs 
(13.6691.4052, 6.6691.1003) and (13.6810.2845, 6.6810.428) represent the same stars 
found twice in overlapping regions of different MACHO fields. The total number of distinct 
stars is therefore 73. We remind the reader that our preferred nomenclature is MACHO* 
followed by the right ascension and declination as listed in Table 1, with no spaces. 

The periods of the RR Lyrae were initially estimated using the IRAF PDM task as 
described in Alcock et al. (1996). The periods, period ratios and modal amplitudes reported 
in Table 1 are the result of a weighted least-squares fit to the RRd lightcurves for the first 
three Fourier terms in each mode and all coupling terms. 

The ratio of first-overtone to fundamental mode amplitudes is shown in Figure 2. It is 
clear that this sample of stars is similar to known RRd's in that, almost without exception. 
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the amplitude of the first-overtone mode is greater than that of the fundamental. We note 
that the selection criteria used above are biased against stars with fundamental periods 
shorter than 0.46 days and that amplitude ratios smaller than unity may therefore be 
under-represented . 

2.4. Absolute Magnitudes from Pulsation Theory 

It is, in theory, possible to estimate absolute magnitudes of RRd stars by combining the 
Stefan-Boltzmann equation and the period-density relations. The result of such an exercise 
is an expression relating period to luminosity, effective temperature, and mass. Unlike 
singly-periodic RR Lyrae stars, the period ratio allows us to estimate the mass based on 
atmospheric pulsation models. Provided that we can determine the effective temperature, 
the set of required observables is complete. Such a relationship is well-known for Cepheid 
variables where the existence of the instability strip boundaries gives rise to an obvious 
period-luminosity relation. For double-mode RR Lyrae stars we expect a similar bounded 
region in the color-magnitude diagram which will also give rise to a period-luminosity 
relationship. 

We will now proceed through the steps necessary to obtain absolute magnitude 
estimates. First, the estimation of masses from periods and period ratios will be described. 
Second, the estimation of effective temperature and the relationship of these stars with 
the fundamental blue edge (FBE) of the instability strip will be described. Third, we 
substitute for both effective temperature and mass in the pulsation equation and produce a 
period-luminosity relation. Finally, we estimate the errors introduced at each step. 
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2.4- 1- Period Ratios and Masses 

We show in Figure 3, the Petersen (1973) diagram of Pi/Po plotted against Pq for 
stars in M15 (Nemec, 1995a), Draco (Nemec, 1995b), IC 4499 (Clement et al, 1986, Walker 
and Nemec, 1996) and the LMC. We additionally mark on the diagram the locus of period 
ratios for model pulsators computed by Cox (1991) using 0PAL91 opacities. 

The distribution of the stars in Figure 3 shows that the LMC stars include a low mass 
extension to the multimode stars discovered in other systems. From Figure 3 it can be 
shown that a relation of the form 

logPo = 0.4521ogA4- 0.228 (1) 

applies across the indicated mass range of the multimode stars in the LMC and the Galactic 
Oosterhoff groups I and II. 

The shortest period of Bailey type a,b stars among the total sample of LMC RR Lyrae 
stars is 0.459 days and, additionally, with a period ratio of about 0.7421 the overtone 
periods are close to one-third of a day and are thus more difficult to discover due to aliasing. 
We mark these period regions as hatched areas in Figure 3 to emphasize the degree to which 
the period-distribution of the multimode stars is affected by selection biases in this work. 
Nevertheless, it is clear that the tight locus of observed period and period-ratio in Figure 3 
for RRd stars born in a wide range of environments is consistent with a tight relationship 
between luminosity, effective temperature and mass. In the following discussion we shall 
assign masses to the LMC stars on the basis of their position in the Petersen diagram 
relative to the Cox (1991) models. We note that there are potential uncertainties in these 
masses due to the lack of specific data on composition of the LMC stars. However, the 
abundance distribution of the LMC field RR Lyrae stars as derived from spectra does not 
deviate strongly from that of the galactic halo ( Alcock et al. 1996). 
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2.4-2. Pulsation Models 

The pulsation equation for model envelopes of RR Lyrae stars has been characterized 
using OPAL95 opacities by Bono et al. (1996) as 

logPo.tr = 11.627 + 0.823 log - 3.506 logTe,t^ - 0.5821ogA1t^ (2) 

over the mass range 0.65-0.75 M0. We have annotated the variables to reflect the assumption 
that the multimode stars are in the transition zone (the region between the fundamental 
blue edge and the first-overtone red edge) of the instability strip. Their equation 2 is a 
linear fit to the results of period determinations for RR Lyrae star models using the most 
recent composition parameters in a study of the non-linear pulsation properties. The fit to 
the pulsation equation is in the spirit of the influential study of RR Lyrae star pulsation by 
van Albada and Baker (1971) and we note in passing that use of their form of the pulsation 
equation does not affect the results and conclusions drawn in this paper. 

We make the further assumption (see, for example, Caputo 1990, Bono and 
SteUingwerf 1994) that there is close similarity between the temperature of the blue 
edge, (FBE) of the fundamental instabihty strip and the transition zone occupied by the 
multimode stars. Ptr is thus taken as the period of the fundamental mode at the blue edge 
of the instability strip. In the range of Ptr of the beat RR Lyrae stars, Tg varies from 
6900 K to 7060 K in excellent agreement with the calculations of the fundamental and 
first-overtone instability strip boundaries (and the zone where both fundamental and first 
overtone are unstable), by Bono and SteUingwerf (1994). Their calculations show that the 
temperature of the strip in which pulsation in either the fundamental or first-overtone mode 
is allowed is relatively insensitive to mass or luminosity and is centered near 6840 K for the 
mass range 0.58-0.75 Mq and luminosity range 1.6-1.9 logi^Q. The assumption that the 
double-mode stars lie near this temperature is also supported by the close similarity of the 
periods of the shortest period LMC Bailey type a,b stars at Pq — 0.457 days to the average 
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period of the multimode stars at Pq — 0.48 days. 

For the FBE, we can derive a relationship between the fundamental period and the 
effective temperature by combining equation 8 of Sandage (1993a) - the relation between 
the shortest period for RRab stars and [Fe/H] - and equation 5 of Sandage (1993b) - the 
relation between logTg and [Fe/H] for RRab stars - to obtain: 

log Te,t^ = 3.816 - 0.0984 log Pt^. (3) 

We see that the FBE is very nearly vertical in the logL, logTe plane. Sandage expresses 
both the FBE period and temperature in terms of abundance [Fe/H] but since the range of 
[Fe/H] in the LMC field RR Lyrae stars is the same as in the calibrating galactic RR Lyrae 
clusters and field, we have written the FBE temperature directly in terms of period. The 
luminosity of the multimode stars is then 

log Ltr = 2.483 + 2.360 log Ptr (4) 

Over the range of the transition periods covering the majority of the LMC multimode stars 
shown in Figure 3, at Pq — 0.46 days, logL^^ = 1.69 and at Pq = 0.52 days, logL^^ = 1.81, 
corresponding to absolute visual magnitudes of My — +0.49 and +0.19 respectively when 
the small bolometric corrections of Kurucz (1979) arc applied. For the longer fundamental 
periods of the multimode stars in the metal-weak galactic globular cluster M 15 with a 
characteristic Pq = 0.54 days, the derived luminosity corresponds to My — +0.05 mag. 
Most other calibrations of RR Lyrae star luminosities have been expressed in terms of 
composition and there has been considerable debate about the slope of the coefficient of 
the term in [Fe/H]. We note that if we adopt the composition value of [Fe/H] = -2.12 for 
M15 ( Buonanno et al. 1989), our cafibration has a roughly similar dependence on [Fe/H] 
to that discussed by Sandage (1993a,b,c) but caution that direct comparison of these 
dependencies is limited by the small number of clusters that contain multimode stars. We 
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also note that our calibration depends on the theoretical pulsation equation which, as we 
will show in the next section, appears to be supported by the luminosities of the LMC 
RR Lyrae stars. Ours is a brighter calibration of RR Lyrae star luminosities than many 
previously derived. This calibration receives support from the independent study of the 
Bailey type c RR Lyrae stars by Simon and Clement (1993) who compared the Fourier 
phase parameters and periods of the first overtone pulsators with linear and hydrodynamic 
pulsation models to determine temperatures and a luminosity calibration. The result of 
their calibration apphed to the Reticulum cluster of the LMC yields roughly the same 
distance modulus as we derive below for the field LMC stars. Their discussion of the galactic 
multimode stars as checks on their calibration of the RRc variables additionally supports 
the assumptions we made above concerning the location of the multimode instability strip 
in the luminosity-temperature plane. 

In his recent monograph. Smith (1995) reviews at least 6 methods, both fully and 
semi-empirical, which yield luminosities roughly 0.4 mags fainter than those derived here. 
Conversely, all these methods would predict the LMC RR Lyrae stars to be too faint by the 
same amount. In this paper we adopt the view that the characteristics of RR Lyrae stars 
in the LMC and the Galaxy are properly comparable and explore the consequences of that 
comparison. 

2.4-3. Error Analysis 

The application of the luminosity calibration derived in the previous section to the 
photometry of our multimode RR Lyrae stars in the LMC necessitates a careful examination 
of the uncertainties. We adopt a mass uncertainty of ±0.05 M.q ( Kovacs et al. 1992) in 
our logPo, logA^ relation. Assuming double-mode RR Lyrae stars lie at the FBE, the 
uncertainty in equation (3), our logTe^j^, log Ptr relation derived from Sandage (1993a,b), is 
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dominated by the uncertainty of the (B-V), logT^^tr relation used to define the temperatures 
of galactic RR Lyrae along the FBE. We return to the unpubhshed model atmosphere 
calculations of Bell as seen in Butler et al. (1978), the model atmospheres of Kurucz 
(1979), and the instability strip analysis of Bono et al. (1996), and by intercomparison 
adopt at the FBE an uncertainty of 0.01 in logTe^^^- Propagating the errors through to 
equation (4), we find that at a given period the uncertainty in logLj^ is ±0.065 resulting 
in an uncertainty in My of ±0.16 mag. We have assumed equation (4) is an exact 
representation and there is negligible error in these bolometric corrections. 

244. The LMC Distance Modulus 

In Figure 1 we show the multimode stars in the V,(V-R) color- magnitude diagram 
(CMD). We adopt the reddening vector with Ry = Av/E(V - R) = 5.0 ( Bcsscll, 1996) 
which is illustrated in Figure 1. Also shown is the adopted photometric uncertainty for each 
point. The insert shows the multimode stars in relation to a representative MACHO LMC 
CMD. Using the temperatures from equation (3), for stars with log^f = 2.6, we adopt the 
(B-V), (V-R) relation ( Bessell 1996), 

{y -R)= 0.004 ± 0.566(5 - V) (5) 

with an uncertainty of ± 0.017 in the zero point. From Butler et al. (1978) we interpolate 
the T (our adopted temperature for the FBE), (B-V) relations between [Fe/H] = -1 and -2 
for [Fe/H] — -1.7, a value appropriate for the LMC RR Lyrae ( Alcock et al. 1996) and 
derive a (V — R)q, logPt,. relation, 

(y -i?)o = 0.19 ±0.15 log Pt, (6) 

with an uncertainty in the zero-point of ±0.02 mag derived from the uncertainties in 
equations (3) and (5) as discussed previously. Adopting our FBE temperature, the intrinsic 
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color, {V — R)o, ranges from 0.14 to 0.15 mag for our multimode stars. Furthermore, we 
note that the mean reddening, < E{V - R) > = 0.049 mag, or < E{B -V) > = 0.086 mag, 
agrees well with other measurements of the reddening toward the LMC ( Bessell 1991). In 
Figure 4 we show the de- reddened V magnitudes plotted against Pq- The error bars noted 
are derived from the uncertainty in V and in (V-R) scaled by the ratio of total to selective 
absorption (Ry = 5), equivalent to ±0.27 mag. The contribution of the uncertainty in our 
adopted intrinsic color is systematic and we have not included it in the uncertainty of each 
individual de-reddened magnitude. Our multimode RR Lyrae have a mean period < Pq > 
= 0.482 corresponding to < My > = 0.37 ± 0.16 mag. The mean de- reddened magnitude 
of our multimode RR Lyrae sample is < Vq > = 18.94 ± 0.03 mag, in excellent agreement 
with the mean de-reddened V magnitude of nearly 180 RR Lyrae found in LMC clusters ( 
Walker 1992). The distance modulus of the LMC naturally follows: (m - M)lmc = 18.57 
± 0.19 mag, where the final error is a combination of the uncertainty in {V — R)o (scaled 
by Rv — 5), the uncertainty in < Vq >, and the uncertainty in < My >■ Fitting the data 
of Figure 4 to the My, logPt^ cahbration with a fixed slope yields the same LMC distance 
modulus. The best fit distance modulus is shown as a horizontal fine in Figure 4. The 
corresponding My, log Pe,tr, calibration is also shown. The distance modulus for the LMC 
multimode RR Lyrae stars derived here is found to be consistent with that (18.5 ± 0.1 mag) 
derived from the galactic calibration of classical Cepheids discussed by Walker (1992). 



3. Discussion 

Walker (1989,1992) has emphasized the apparent discrepancy of 0.3 magnitudes 
between RR Lyrae star luminosities derived from Baade-Wesselink analyses of galactic stars 
and the luminosities of LMC cluster variables if the Cepheid distance calibration is correct. 
We have shown here that there is a consistency between the luminosities derived from 
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the pulsation equation, and independent luminosity measures such as the LMC Cepheid 
distance cahbration and the SN1987A ring distance ( Gould 1995). Consequently, the 
pulsation properties of the multimode RR Lyrae stars appear to be a useful tool in distance 
studies through their calibration of the more common and easily discovered RRab stars. 

The unusually large range in Pi/ Pq and the range in implied masses appear to present 
new questions about the evolution from the ZAHB. The Lee and Demarque (1990) and 
Sweigert (1987) HB evolutionary tracks do not carry the evolution of stars with masses as 
low as O.55M0 as far as the instability strip but it does not appear that such stars could 
be of similar luminosity to the HB when they reach the instability strip following redward 
evolution. Cox (1991) goes further and states that stars with mass of O.SSM© would 
never reach the instability strip at all. However as we see in Figure 3, there do exist LMC 
multimode RR Lyrae stars with this mass as derived from pulsation theory. In Alcock et 
al. (1996), we argued that the field HB of the LMC has a red morphology and that as a 
result, the beat RR Lyrae stars are rare. Our more recent detections still show fewer of 
these stars among the LMC field population than, say in the Ooll cluster M15. We remind 
the reader that our survey is incomplete and that this fraction may rise. In the same paper, 
we also argued that the LMC field RR Lyrae star abundance distribution was similar to the 
halo galactic distribution with a modal [Fe/H] = —1.7. These findings led us to conclude 
that the LMC field was younger than most galactic globular clusters. 

An implication of Figure 4 and the period-luminosity relation contained in equation (4) 
is a differential rescaling of the distances to galactic globular clusters (and the galactocentric 
distance) when HB luminosities are the criteria. Oosterhoff group II clusters are the most 
affected. 

Walker (1992), in an important paper, supposed that the distance to the LMC was 
correctly defined by the galactic Cepheid distance scale and inferred that the LMC RR 
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Lyrae stars were 0.3 mag brighter than previously suggested. He examined the consequences 
of a global increase in HB star luminosities on the age and distance of galactic globular 
clusters. In this paper we have determined the RR Lyrae star luminosities to be significantly 
brighter than previously thought and have shown that this revision is most important for 
the OoII clusters in the Galaxy. For illustrative purposes we discuss the archetypal Galactic 
cluster M15 which has multimode RR Lyrae stars with a mean fundamental period of Pq 
— 0.54 days. Equation (4) implies the HB My — +0.05 mag. Buonanno et al. (1989) hst 
the brightness difference between MSTO and HB as 3.54 mag for this cluster, leading to 
Mv,to = 3.59 mag. These authors extract from evolutionary calculations the relation 

logtg = -0.41 + 0.37My,to - 0.43y - 0.13[Fe/i/] (7) 

between age, turn-off luminosity and composition. For M15 we take Y — 0.23 and [Fe/H] 
= -2.15 (see Buonanno et al. 1989)). This equation yields an evolutionary age of 12.6 
Gyr. (This estimate does not include the effect of any alpha-particle enrichment.) Minniti 
et al. (1996) have found [0/Fe] = 0.45 for M15 and we estimate that this results in a 
reduction in age of 0.6 Gyr based on the work of Safaris et al. (1993). The age estimate 
of 12.0 Gyr has the characteristic uncertainty of ±1.5 Gyr associated with MSTO ages. 
These age estimates can undoubtedly be refined by the fitting of evolutionary tracks to the 
recalibrated CMD's. We note that compared to a HB absolute magnitude of My — -1-0.49 
mag for M15 ( Armandroff 1989) our calibration produces a reduction in age of 33% which 
is relatively independent of the assumed value of alpha-element enrichment. 

A fuller exploration of the parameter space of cosmological models defined by the age 
of the oldest stars is beyond the scope of this paper but we note that a 33% reduction in 
the estimated age of the oldest stars will help reconcile globular cluster ages with recent 
estimates of the Hubble constant, Hq. 
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Fig. 1. — Observed median V magnitudes are plotted against derived (V-R) colors for the 
multimode RR Lyrae stars. The reddening vector is marked with an arrow. The insert shows 
the CMD of the multimode stars (filled circles), in relation to a representative CMD in an 
LMC field. 
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Fig. 2. — The ratio of the amphtudes of the first Fourier components for the first-overtone and 
fundamental modes, Ai/Aq, are plotted against the fundamental mode period Pq. Almost 
without exception, the amplitude of the first overtone mode is greater than that of the 
fundamental - a property that this sample shares with previously discovered RRd stars in 
globular clusters and dwarf spheroidal galaxies. We caution, however, that our sample is 
incomplete due to period-dependent selection biases and that more RRd with Ai/Aq < 1 
may yet be discovered in our LMC sample. 
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Fig. 3. — The Petersen diagram for multimode RR Lyrae stars where the first-overtone to 
fundamental period ratios Pi/Pq are plotted against Pq. As described in the text, the LMC 
stars all have Pq < 0.5 days while the longer period stars are from M15 and Draco. The loci 
of model pulsators calculated by Cox (1991) are shown for the masses of 0.55, 0.65 and 0.75 
M.Q. The near- vertical line at Pq — 0.46 days marks the location of the one-third day alias 
for the first-overtone period. 
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Fig. 4. — The absorption-corrected magnitudes, Vq, are plotted against the fundamental 
period Pq- The hne has the slope given by the equation (4). 



